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Abstract 

Beginning with a Lagrangian that is consistent with both Newtonian gravity and the momentum-velocity 
relation of special relativity, an approximate relativistic orbit equation is derived that describes relativistic 
corrections to Keplerian orbits. Specifically, corrections to a Keplerian orbit due to special relativity include: 
precession of perihelion, reduced radius of circular orbit, and increased eccentricity. The prediction for the 
rate of precession of perihelion of Mercury is in agreement with existing calculations using only special 
relativity, and is one sixth that derived from general relativity. All three of these corrections are qualitatively 
correct, though suppressed when compared to the more accurate general-relativistic corrections in this limit. 
The resulting orbit equation has the same form as that derived from general relativity and is easily compared 
to that describing Kepler's orbits. This treatment of the relativistic central-mass problem is complementary 
to other solutions to the relativistic Kepler problem, and is approachable by undergraduate physics majors 
whom have not had a course dedicated to relativity. 

PACS numbers: 04.25.Nx, 95.30.Sf, 03.30.+p, 95.10.Ce 
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I. INTRODUCTION 



The relativistic contribution to the rate of precession of periheUon of Mercury is calculated 
accurately using general relativity [1-6]. However, the problem is commonly discussed in un- 
dergraduate and graduate classical mechanics textbooks, without introduction of an entirely new, 
metric theory of gravity. One approach [7-11] is to define a Lagrangian that is consistent with 
both Newtonian gravity and the momentum- velocity relation of special relativity. The resulting 
equation of motion is solved perturbatively, and an approximate rate of precession of perihelion of 
Mercury is extracted. This approach is satisfying in that all steps are proved, and a brief introduc- 
tion to special relativity is included. On the other hand, one must be content with an approximate 
rate of precession that is about one-sixth the correct value. Another approach [12-16] is that of a 
mathematical exercise and history lesson. A modification to Newtonian gravity is given, without 
proof, resulting in an equation of motion that is the same as that derived from general relativity. 
The equation of motion is then solved using appropriate approximations, and the correct rate of 
precession of perihelion of Mercury is extracted. Both approaches provide an opportunity for stu- 
dents of classical mechanics to understand that relativity is responsible for a small contribution to 
perihelic precession and to calculate that contribution. 

We present a review of the approach using only special relativity, followed by an alternative 
solution of the equation of motion derived from Lagrange's equations. An approximate rate of 
perihelic precession is derived that agrees with established calculations. This effect arises as one 
of several small corrections to Kepler's orbits, including reduced radius of circular orbit and in- 
creased eccentricity. The method of solution makes use of coordinate transformations and the 
correspondence principle, rather than the standard perturbative techniques, and is approachable by 
undergraduate physics majors. 

A relativistic particle of mass m orbiting a central mass M is commonly described by the 
Lagrangian [7-11, 17-22] 

L = -mc^-f-^ - U{r), (1) 

where: 7"^ = -^1 — v^/c^; v'^ = f'^ + and U{r) = —GMm/r. (G is Newton's universal 
gravitational constant, and c is the speed of light in vacuum.) The equations of motion follow from 
Lagrange's equations, 

ddL dL „ 



for each of {qi} — {9, r}, where (ji = dqi/dt. The results are 

d 



^^hr'e]^0, (3) 



which implies that i = 'yr 9 = constant; and 



7r + 7r + - -fr9^ = 0. (4) 



The first of these [Eq. (3)] is the relativistic analogue to the Newtonian equation for the conserva- 
tion of angular momentum per unit mass, and is used to eliminate 9 in Eq. (4), 

ir9' (5) 

Time is eliminated by successive applications of the chain rule, together with the conserved angular 
momentum [23-27]: 

r = — Tq-; (6) 
^ a9 r 

and, therefore, 

7r = -7r 27^~- 

7r"^ dt/^ r 

Substituting Eqs. (5) and (7) into the equation of motion Eq. (4) results in 

1 

^"^2--GMi + -=^. (8) 
a.9'^ r r 

We anticipate a solution of Eq. (8) that is near Keplerian and introduce the radius of a circular 
orbit for a nonrelativistic particle with the same angular momentum, = i'^/ GM. The result is 

''''■"+'-^ = l + A, (9) 



^Q2 J, J. 

where A = 7 — 1 is a velocity-dependent correction to Newtonian orbits due to special relativity. 
The conic-sections of Newtonian mechanics [28-33] are recovered by setting A = (c — > 00): 

d^ Tc , Tc 



d6'2 r r 

which implies that 



+ - = 1, (10) 



- = l + ecos^, (11) 
r 



where e is the eccentricity. 
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n. KEPLERIAN LIMIT AND ORBIT EQUATION 



The planets of our solar system are described by near-circular orbits (e <^ 1) and require 
only small relativistic corrections {v/c <^ 1). Mercury has the largest eccentricity (e ^ 0.2), 
and the next largest is that of Mars (e fa 0.09). Therefore, A [defined after Eq. (9)] is taken 
to be a small relativistic correction to near-circular orbits of Newtonian mechanics (Keplerian 
orbits). Neglecting the radial component of the velocity, which is small compared to the tangential 
component, and expanding 7 to first order in v^/c^ results in 

X^\{r9/cf. (12) 

(See Sec. IV and App. A for a thorough discussion of this approximation.) Once again using the 
angular momentum to eliminate ^, this may be expressed as A ^ |(£/rc)^7~^ fti |(£/rc)^(l — 2A), 
or 

X^\{^/rcf. (13) 
The equation of motion Eq. (9) is now expressed approximately as 

1 /rA^ 

H fa 1 + 



d9^ r r 2 \r 

where e = {GM/lcf. The conic-sections of Newtonian mechanics, Eq. (10) and Eq. (11), are 
now recovered by setting e = (c — > 00). The solution of Eq. (14) for e 7^ approximately 
describes Keplerian orbits with small corrections due to special relativity. 

If e is taken to be a small relativistic correction to Keplerian orbits, it is convenient to make the 
change of variable 1/s = rg/r — 1 -C 1. The last term on the right-hand-side of Eq. (14) is then 
approximated as (rc/r)^ ^l-\-2/s, resulting in a linear differential equation for l/s{9): 

2 1 2(1 -e)l ^ 
e dc^^ s e s 

The additional change of variable a = 6y/l — t results in the familiar form: 

H fa 1, (16) 



dct^ s s 

where Sc = 2(1 — e)/e. The solution is similar to that of Eq. (10) 



I + Acqs a, (17) 
s 
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where A is an arbitrary constant of integration. In terms of the original coordinates, Eq. (17) 
becomes 



where 



— 1 + 6 008^6', (18) 

r 



1 - e 

Tc = To- ^, (19) 

hA 

e=-^, (20) 
1 



K=(l-e)2. (21) 

According to the correspondence principle, Kepler's orbits, Eq. (11), must be recovered in the 
Umit e — > (c — > oo), so that ^€A = eis the eccentricity of Newtonian mechanics. To first order 
in e, Eqs. (19)-(21) become 

fe«ire(l-|e), (22) 
e«e(l + |e), (23) 
K ^ 1 - |e, (24) 

so that relativistic orbits in this limit are described concisely by 

''"^^'^^^ ^l + e{l + |e) cos (1 - le)0. (25) 

This approximate orbit equation has the same form as that derived from general relativity [34], 



"-^ — ^ ^l + e{l + 3e) cos (1 - 3e)e, 



(26) 

r 

and clearly displays three characteristics: precession of perihelion; reduced radius of circular orbit; 
and increased eccentricity. 



III. CHARACTERISTICS OF NEAR-KEPLERIAN ORBITS 

The approximate orbit equation Eq. (25) predicts a shift in the perihelion through an angle 

= 27r{k-^ -l)^7re (27) 

per revolution. This prediction is identical to that derived using the standard approach [7-11, 35] to 
incorporating special relativity into the Kepler problem, and is compared to observations assuming 
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that the relativistic and Keplerian angular momenta are approximately equal. For a Keplerian orbit 
[28-33] e = GMa{l - e^), where G = 6.670 x IQ-^^m^kg-^ s'^, M = 1.989 x 10=^° kg is the 
mass of the Sun, and a and e are the semi-major axis and eccentricity of the orbit, respectively. 

Therefore, the relativistic correction defined after Eq. (14), 

GM 



(28) 



is largest for planets closest to the Sun and for planets with very eccentric orbits. For Mercury 
[36, 37] a = 5.79 x 10^° m and e = 0.2056, so that e 2.66 x 10"^ (The speed of light is taken 

to be = 8.987554 x 10^^ m^/s^.) According to Eq. (27), Mercury precesses through an angle 

~ o ON = 8-36 X 10-^ rad (29) 

&a{l — e^j 

per revolution. This angle is very small and is usually expressed cumulatively in arc seconds per 
century. The orbital period of Mercury is 0.24085 terrestrial years, so that 

100 yr 360 x 60 x 60 , ^ , 

AO = — — X X (30a) 

0.24085 yr 27r ^ ' 

ft! 7.16arcsec/century. (30b) 

The general relativistic (GR) treatment results in a prediction of 43.0 arc sec/century [34—55], 
and agrees with the observed precession of perihelia of the inner planets [37^6, 56-59]. His- 
torically, this contribution to the precession of perihelion of Mercury's orbit precisely accounted 
for the observed discrepancy, serving as the first triumph of the general theory of relativity [1^]. 
The present approach, using only special relativity, accounts for about one-sixth of the observed 
discrepancy, Eq. (30b). Precession due to special relativity is illustrated in Fig. 1. 

The approximate relativistic orbit equation Eq. (25) [or Eq. (18) together with Eqs. (22)-(24)] 
predicts that a relativistic orbit in this limit has a reduced radius of circular orbit (e = 0). This 
characteristic is not discussed in the standard approach to incorporating special relativity into the 
Kepler problem, but is consistent with the GR description. An effective potential naturally arises 
in the GR treatment of the central-mass problem [34, 36, 40, 41, 43, 46, 53, 54, 60, 61], 

GM e GMf 

Kff = + - (31) 

r Zr c r 

that reduces to the Newtonian effective potential in the limit c — > oo. In the Keplerian Umit, 
the GR angular momentum per unit mass I is also taken to be approximately equal to that for a 
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Keplerian orbit [34, 3>6-A3, 46, 54, 61, 62]. Minimizing V^ff with respect to r results in the radius 
of a stable circular orbit, 

Rc = \n + ^rcVl - 12e « rc(l - 3e), (32) 

so that the radius of circular orbit is predicted to be reduced: Rc — ~ — 3erc. (There is also 
an unstable circular orbit. See Fig. 2.) This reduction in radius of circular orbit is six times that 
predicted by the present treatment using only special relativity, for which fc — rc ~ — |erc. [See 
Eq. (22).] 

Most discussions of the GR effective potential Eq. (31) emphasize relativistic capture, rather 
than reduced radius of circular orbit. The 1/r^ term in Eq. (31) contributes negatively to the 
effective potential, resulting in a finite centrifugal barrier and affecting orbits very near the central 
mass (large-velocity orbits). (See Fig. 2.) This purely GR effect is not expected to be described by 
the approximate orbit equation Eq. (25), which is derived using only special relativity and assumes 
orbits very far from the central mass (small-velocity orbits). 

An additional characteristic of relativistic orbits is that of increased eccentricity. Equation (25) 
predicts that a relativistic orbit will have increased eccentricity, when compared to a Keplerian 
orbit with the same angular momentum: e — e fa |ee. [See Eq. 23.] This characteristic of 
relativistic orbits also is not discussed in the standard approach to incorporating special relativity 
into the Kepler problem, but is consistent with the GR description. The GR orbit equation in this 
Keplerian limit Eq. (26) predicts an increase in eccentricity e — e ^ See, which is six times that 
predicted by the present treatment using only special relativity. 



IV. DISCUSSION 



The approximate orbit equation in Eq. (25) provides small corrections to Kepler's orbits due 
to special relativity. [Compare Eqs. (25) and (11).] A systematic verification may be carried out 
by substituting Eq. (25) into Eq. (14), keeping terms of orders e, e, and ee only. The domain of 
validity is expressed by subjecting the solution Eq. (25) to the condition 

- - K 1 (33) 
r 

for the smallest value of r. Evaluating the orbit equation Eq. (25) at the perihelion (rp) results in 

Tn 1- ke 
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The substitution of this result into Eq. (33) results in the domain of validity: 

e(l + |e) + e < 1. (35) 

Therefore, the relativistic eccentricity e — e{l + |e) <^ 1, and Eq. (25) is limited to describing 
relativistic corrections to near-circular (Keplerian) orbits. Also, the relativistic correction e <^ 1, 
and thus the orbit equation Eq. (25) is valid only for small relativistic corrections. 

In Sec. II the relativistic correction to Keplerian orbits A [defined after Eq. (9)] is approximated 
by: neglecting the radial component of the velocity, = + (r^)^ ^ (r^)^; and keeping 
terms only up to first order in the expansion ^ 1 + ^{v/cy. Neglecting the radial component 
of the velocity in the relativistic correction A is consistent with the assumption of near-circular, 
approximately Keplerian orbits, and is complementary to the assumption rc/r — 1 <^ 1 preceding 
Eq. (15). It is emphasized that the radial component of the velocity is neglected only in the 
relativistic correction A; it is not neglected in the derivation of the relativistic equation of motion 
Eq. (9). That there is no explicit appearance of r in the relativistic equation of motion Eq. (9), 
other than in the definition of 7, is due to a fortunate cancellation after Eq. (7). Furthermore, the 
approximate orbit equation Eq. (25) has the same form as that [Eq. (26)] arising from the GR 
description, in which the radial component of the velocity is not explicitly neglected. 

A first-order series approximation for 7 is consistent with the Keplerian limit. Interestingly, 
however, the problem is soluble without truncating the series. (See Appendix A.) This more 
elaborate derivation yields diminishing returns due to the complementary assumption Vc/r — 1 ^ 
1, which implicitly constrains the velocity to be much smaller than the speed of light. The resulting 
approximate orbit equation [Eq. (A13)], 

r (1- -e) 

1 + e(l + f e) cos (1 - \e)9, (36) 

is almost identical to that derived using a much simpler approach in Sec. II [Eq. (25)]. Although 
this alternative orbit equation predicts a relativistic correction to eccentricity that is nearly one- 
half that of the GR result, it lacks the symmetry of the GR orbit equation Eq. (26), and it does not 
provide significant further qualitative understanding of relativistic corrections to Keplerian orbits. 

V. CONCLUSION 

The present approach to incorporating special relativity into the Kepler problem results in an 
approximate orbit equation [Eq. (25)] that has the same form as that derived from general rela- 
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tivity in this limit [Eq. (26)] and is easily compared to that describing Kepler's orbits [Eq. (10)]. 
This orbit equation clearly describes three corrections to a Keplerian orbit due to special relativ- 
ity: precession of perihelion; reduced radius of circular orbit; and increased eccentricity. The 
predicted rate of precession of perihelion of Mercury is identical to established calculations us- 
ing only special relativity. Each of these corrections is exactly one-sixth of the corresponding 
correction described by general relativity in the Keplerian limit. 

This derivation of an approximate orbit equation is complementary to existing calculations of 
the rate of precession of perihelion of Mercury using only special relativity. The central-mass 
problem is described by a Lagrangian that is consistent with both Newtonian gravity and the 
momentum- velocity relation of special relativity. However, coordinate transformations and the 
correspondence principle are used to solve the equations of motion resulting from Lagrange's 
equations, rather than the standard perturbative methods. The resulting closed-form, approximate 
orbit equation exhibits several characteristics of relativistic orbits at once, but is limited to de- 
scribing small relativistic corrections to approximately Newtonian, near-circular orbits. This orbit 
equation, derived using only special relativity, provides a qualitative description of corrections to 
Keplerian orbits due to general relativity. Exact solutions to the special relativistic Kepler prob- 
lem require a thorough understanding of special relativistic mechanics [63, 64] and are, therefore, 
inaccessible to most undergraduate physics majors. The present approach and method of solution 
is understandable to nonspecialists, including undergraduate physics majors whom have not had a 
course dedicated to relativity. 



Appendix A 

For near-circular orbits the radial component of the velocity is small compared to the tangential 
component and is neglected in the relativistic correction A = 7 — 1. [See Eq. (9).] An exact 
infinite-series representation, 

^ (2ii)! (reV' 

is used, rather than the approximate first-order truncated series used in Sec. II. (The remainder of 
this derivation parallels that of Sec. II.) Conservation of angular momentum is used to eliminate 9, 
resulting in 




(A2) 
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in terms of which the equation of motion Eq. (9) is now expressed approximately as 

'^^ + !:£~1 + A. (A3) 



The conic-sections of Newtonian mechanics, Eqs. (10) and (11), are recovered by setting c — )■ oo. 

For near-circular, approximately Keplerian orbits it is convenient to make the change of variable 
1/s = Tc/r — 1 <^ 1, so that (rc/r)^° 1 + 2n/s. In terms of this new variable, the relativistic 
correction Eq. (A2) is 

\^{^^-l) + e^l-, (A4) 
s 

where 7^ = 1/ \/l — e, and the following series is identified: 



2n(2n)! / c x 3 
^2^1,^) ^^^^ 



2n 



Thus, the equation of motion Eq. (A3) is linearized: 



+ ^- ~ 1. (A6) 



7g — 1 d9^ s 7g — 1 s 
The additional change of variable a = 6' i^/T^^^Tf results in the familiar form: 

d^ Sr s, 



+ - fti 1, (A7) 



do;^ s s 

where Sc= (1 — ^ll)/{le ~ !)• Th^ solution is similar to that of Eq. (10) 



s 



1 + A COS a, (A8) 



where A is an arbitrary constant of integration. In terms of the original coordinates, Eq. (A8) 
becomes 

— f^l + ecosK^, (A9) 

r 

where (including first-order approximations) 

fc = ^c^^^3^^ ^ n{l - \e), (AlO) 

e = - + -A (All) 

k={l-efS^ ^l-\e. (A12) 
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According to the correspondence principle, Kepler's orbits [Eq. (11) with e < 1] must be recovered 
in the limit e — > (c — > oo), so that |eA = e is the eccentricity of Newtonian mechanics. 
Therefore, relativistic orbits in this limit are described concisely by 

r (I — ^e) 

^ ^ ^ ?^l + e(l + fe)cos(l-|e)g. (A13) 

This alternative orbit equation differs from that [Eq. (25)] derived using the much simpler approach 
in Sec. II only in the relativistic correction to eccentricity. See the discussion at the end of Sec. IV. 
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FIG. 1: Relativistic orbit in a Keplerian limit (solid line), as described by Eq. (25), compared to a corre- 
sponding Keplerian orbit (dashed line) [Eq. (1 1)]. The precession of perihelion is one orbital characteristic 
due to special relativity and is illustrated here for < < 207r. This characteristic is exaggerated by 
both the choice of eccentricity (e = 0.3) and relativistic correction parameter (e = 0.2) for purposes of 
illustration. Precession is present for smaller (non-zero) reasonably chosen values of e and e as well. (The 
same value of e is chosen for both curves.) 
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FIG. 2: Effective potential commonly defined in the Newtonian limit to General Relativity (solid line) 
[Eq. (31)], compared to that derived from Newtonian mechanics (dashed line). The vertical dotted lines 
identify the radii of circular orbits, Rc and r^, as calculated using general relativity and Newtonian me- 
chanics, respectively. General relativity predicts [Eq. (32)] a smaller radius of the circular orbit than that 
predicted by Newtonian mechanics. The value e = 0.06 is chosen for purposes of illustration. 
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